Strong Coupling to Two-Dimensional Anderson Localized Modes 
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We use a scattering formalism to derive a condition of strong coupling between a resonant scatterer 
and an Anderson localized mode for electromagnetic waves in two dimensions. The strong coupling 
regime is demonstrated based on exact numerical simulations, in perfect agreement with theory. The 
strong coupling threshold can be expressed in terms of the Thouless conductance and the Purcell 
factor, thus connecting key concepts in transport theory and cavity quantum electrodynamics. 
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Enhancing and controlling light-matter interaction has 
been an issue of tremendous interest for years. The pio- 
neering prediction of the dependance of the spontaneous 
decay rate of an emitter on its environment in the weak 
coupling regime, known as the Purcell effect [l| , was ob- 
served in optics by Drexhage The development of 
cavity quantum electrodynamics (QED) has led to the 
observation of the strong coupling regime, characterized 
by Rabi oscillations of the excited-state population, or 
a splitting in the frequency spectrum Strong cou- 
pling has been demonstrated with single atoms in en- 
gineered vacuum cavities and in condensed matter 
using quantum-well or quantum-dot excitons in micro- 
cavities or photonic crystals [H-Q. In nanophotonics, 
surface-plasmon modes on metallic nanoparticles or sub- 
strates provide subwavelength light confinement without 
a physical cavity, and strong coupling has been reported 
with quantum dots or molecules [8l4l2l|. Multiple scat- 
tering in disordered media provides an alternative route 
since confined modes can be produced by the mecha- 
nism of Anderson localization [Tsj. Large Purcell fac- 
tors have been measured using quantum dots and An- 
derson localized modes in disordered photonics crystal 
waveguides [l^ . and the strong coupling regime is ex- 
pected to be within experimental reach [16| . In these one- 
dimensional structures, even fabrication imperfections in 
otherwise perfect waveguides generate efficient localiza- 
tion on the micrometer scale [15|. In addition to multi- 
ple scattering, near-field interactions in disordered media 
also contribute to an enhancement of light-matter inter- 
action at scales below the wavelength, as demonstrated 
in the weak-coupling regime [iTj . 

In this Letter, we study the strong coupling regime be- 
tween a resonant dipole scatterer and a two-dimensional 
(2D) Anderson localized mode. Based on a scattering 
formalism for electromagnetic waves, we use numerical 
simulations to demonstrate the strong coupling regime in 
2D localized systems, in perfect agreement with a simple 
mode coupling theory. We show that the strong cou- 
pling criterion can be expressed in terms of the Thouless 
conductance and the Purcell factor. This establishes an 



interesting connection between concepts in transport the- 
ory and cavity QED, that might be useful for the design 
and/or analysis of future experiments. 

In a first part, we use the LDOS spectrum to charac- 
terize an Anderson localized mode. We consider a two- 
dimensional disordered medium and Transverse Electric 
(TE) waves (electric field perpendicular to the plane con- 
taining the 2D scatterers), so that we are left with a scalar 
problem. To introduce the methodology, let us first con- 
sider the canonical situation of a non-absorbing environ- 
ment placed in a closed cavity. In this case, one can 
define an orthonormal discrete basis of eigenmodes with 
eigenfrequencies w„ and eigenvectors e„(r). The electro- 
magnetic response of the medium can be expanded over 
the set of eigenmodes [ll| 



G(r,r',c.) = ^, 
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where c is the speed of light in vacuum, uj the fre- 
quency and G(r, r'jLij) the Green function. The latter 
connects a source dipole with dipole moment p located 
at r' to the radiated field at point r through the rela- 
tion E(r) = /ioa;^G(r, r',aj)p. In the general case of a 
leaky system the weak losses out of each mode can be 
taken into account phenomenologically using a linewidth 
r„. The electric part of the LDOS, relevant for the 
coupling with electric dipoles, is defined as p(r, uj) = 
2a;/(7rc2)ImG(r,r,w) 0. Therefore the LDOS spec- 



trum is given by 
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where An — |e„(r)p. The LDOS spectrum contains all 
the relevant parameters of a given mode (central fre- 
quency, linewidth and local intensity). A major inter- 
est is that it can in principle be determined experimen- 
tally from fluorescent lifetime measurements, even at the 
nanoscale in complex geometries [205. For convenience, 
we also define the Purcell factor associated to a given 
mode n, and a position r, as Fp = p„(r, aj,i)/po, where 
po = w/(27rc^) is the vacuum LDOS in 2D. 

In order to investigate Anderson localization numer- 
ically, we consider an assembly of 2D point scatter- 
ers described by their electric polarizability a(uj) = 



2 



(2ro/fcg)(a;o — w — iTo/2)~^, where fcg = w/c, wq is 
the resonance frequency and Tq the hnewidth. This 
form of the polarizabihty satisfies the optical theorem 
(energy conservation). In the present work, we have 
fixed luq — 3 X 10^^ s^^ (visible optical radiation) and 
To = 5 X 10^® s^^ Wo- With such parameters (wide 
resonance) , the scattering cross-section of the scatterers 
is constant over the spectral range we will consider below. 
The scatterers are randomly distributed in a cylinder of 
radius R. In order to compute the LDOS at the central 
point r5, we need to compute the field scattered at 
when the system is illuminated by a source dipole p also 
located at rs [see Fig. [Tfa)]. The exciting field on scat- 
terer number i is given by the self-consistent equation 

^2 

(3) 

where is the position of scatterer number i. The 2D 
vacuum Green function is Go(r,r',w) — (i/4) HQ^''(fco|r — 
r'l), where Hq^^ is the zero-order Hankel function of the 
first kind. For a system with TV scatterers, the linear sys- 
tem of N self-consistent equations can be solved numeri- 
cally. Once the exciting field on each scatterer is known, 
it is possible to compute the scattered field at rs, and to 
deduce the Green function and the LDOS p(r5,a;). 

Let us consider one configuration of the random sys- 
tem, with iV = 5000 scatterers in a cylinder of radius 
i? = 20 /im. Two computed LDOS spectra, with the 
same bandwidth but centered on two different central 
frequencies wj? = 4 x 10^^ s~^ (diffusive regime) and 
= 8.7 X 10^"' s^^ (localized regime), are shown in 
Fig. [T](b) and Hfc), respectively. To choose these two 
frequencies, we have estimated the localization length 
by ^ = exp[7rRe (fccff)^s/2], with tbe scattering 
mean freepath and fceS the effective wavenumber in the 
medium [2l|, . For a rough estimate, we have made the 
approximation k^s ~ k^ + i/ (2£s), valid in the weak scat- 
tering limit. In the spectrum shown in Fig.[IJb), one has 
^ « 2 X lO^i? and the sample is in the diffusive regime. 
We observe a pretty smooth profile corresponding to the 
intuitive picture of a continuum of modes. Conversely, 
in Fig. [TJc), the localization length is ^ « R/38 and the 
sample is in the localized regime. We observe very sharp 
and well-separated peaks, each of them being a signa- 
ture of a localized mode. A zoom on one of the LDOS 
peaks, as displayed in Fig. [T](d), shows that it can be 
perfectly fitted by a Lorentzian lineshape as in Eq. 
demonstrating the relevance of this description. This iso- 
lated Anderson localized mode will be denoted by mode 
M in the following, and will be used to demonstrate nu- 
merically the strong coupling regime. It is characterized 
by an eigenfrequency um « 8.67 x 10^^ s~^, a linewidth 
Tm = 9.24 X lO^s^i and a Purcell factor Fp w 140. 

In the second part, we now describe theoretically the 
coupling between a resonant dipole scatterer and an An- 
derson localized mode, in order to establish the threshold 
condition for strong coupling. In very general terms, the 
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FIG. 1: (Color online) (a) Sketch of the system, (b) LDOS 
spectrum centered at ~ 4 x 10^'' (diffusive regime), (c) 
LDOS spectrum centered at oj^ — 8.7 x 10^* s~^ (local- 
ized regime), (d) Zoom on one peak in the localized regime. 
Circles correspond to a fit by Eq. 0. 

scattering medium is described by the scattered Green 
function S(r, r', w) = G(r, r', w)— Go(r, r', cj), where Gq is 
the vacuum Green function (or more generally the Green 
function in a reference medium). The resonant dipole 
scatterer, placed at position rg, is described by its elec- 
tric polarizabihty as in vacuum (or in the same reference 
medium), such that its induced dipole moment reads 

p(w) = eoas(w)E'='''=(rs,w), (4) 

where E°^'^{rs,uj) is the exciting field. The eigenmodes 
of the coupled systems are found by assuming that the 
exciting field is provided by the polarizable scatterer itself 
(no external illumination), so that [1^ 

E-^^rs,oj) = fioio^S{rs,rs, uj)p{uj). (5) 

Combining Eqs. ^ and (O, one obtains the implicit 
equation satisfied by the eigenfrequencies of the coupled 
system 

^2 

^as(w)S(rs,rs,w) = 1. (6) 

Let us stress that this relation is very general and rules 
the coupling between the scatterer and its environment 
whatever the strength of this coupling. 
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In the case of an Anderson localized mode centered at 
lomi the Green function in the vicinity of ujm is given by 



e*M{r')eM{r) 



2ujM UJM -u} - ITm/^ 



(7) 



We consider a scatterer resonant with mode M, described 
by a polarizability 



asiuj) 



2c' 
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where Tg and r;^^ are, respectively, the radiative and 
intrinsic non-radiative decay rate of the scatterer. Note 
that Tg also appears in the numerator as it includes the 
oscillator strength of the polarizable scatterer. The com- 
plex eigenfrequencies Q of the coupled system are solu- 
tions of Eq. (|6]), in which as{io) is given by Eq. (jS]) and 
S(r5, r5, uj) is deduced from Eq. ([7]). One finds two solu- 
tions 
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where gc = iTgTM^p/^Y''^ is the classical coupling con- 
stant. Under the condition 



> 
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(10) 



the two new eigenmodes of the coupled system are no 
longer degenerated. This defines the strong coupling 
regime. The frequency splitting between the two eigen- 
modes is given by the Rabi frequency fifl, such that 



9c 



16 



1/2 



(11) 



The spectral width F of the new eigenmodes is the av- 
erage of r^-'^ and Fm, i.e., of the intrinsic linewidths of 
the uncoupled system. Let us note that condition ((TU)) is 
not sufficient to ensure that the Rabi splitting is larger 
than the linewidth (this would be a necessary condition 
to observe Rabi oscillations in the time domain). One 
needs to satisfy the more restrictive condition 2VLr > F, 
that reads 



9. 
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The scattering formalism developed here encompasses 
the case of a quantum emitter with transition dipole 
d coupled to a localized mode. Replacing Tg by the 
spontaneous decay rate F^ = ^oujf^(P/{2h) in Eq. ([TU|. 
and assuming F^^ = 0, one recovers the condition 
QmAm/ojm > fi^o/{8(f ) that is given, e.g., in Ref. p^ . 
Qm being the quality factor of the localized mode. Fi- 
nally let us note that the weak-coupling regime is recov- 
ered when Eq. ([TU]) is not satisfied, in the limit Fm ^ F5. 



In this limit, Eq. ([9]) shows that the system remains de- 
generated. The resonance of the scatterer broadens with- 
out affecting the localized mode. The broadening in this 
regime (or the change in the spontaneous decay rate for a 
quantum emitter) is exactly given by the Purcell factor. 

The expected strong coupling regime can be checked 
using exact numerical simulations. We consider the same 
system as in Fig. [IJa) , and add at position rs a resonant 
dipole scatterer (probe scatterer), tuned to the resonance 
frequency ujm of the localized mode M identified in the 
spectrum in Fig. [IJc). The polarizability of the probe 
scatterer is given by Eq. ([H]), with F^^^ = 0. A sketch 
of the system is represented in Fig[2l^a). Under an il- 
lumination by an external field Eo(r, w) (plane wave), a 
system of N self-consistent equations similar to can 
be written 



, ,2 



H — 5-as(w)Go(ri,rs,w)Es, (13) 

where the exciting field E5 on the probe scatterer is given 
by 

2 ^ 

Es = Eo(rs,a;) + ^^a(w) ^ Go(rs, r^, cj)Ej. (14) 

3=1 

Solving this linear systems with A'^ -I- 1 equations al- 
lows us to compute the induced dipole moment of the 
probe scatterer ps(w) = eoas(i^)E5(r5, w). We show 
in Fig. [2jb) the resulting spectrum for five different val- 
ues of the radiative linewidth F^ (increasing from top 
to bottom). The Rabi splitting 2Q,ii increases with F;^, 
as expected from theory since the coupling strength gc 
scales as (Fg)^/^. The dependence of the Rabi splitting 
on Fg extracted from the numerical simulations is shown 
in Fig. [Jfc) . Excellent agreement is found with the the- 
ory. In summary, the simulations, performed without 
any approximation, have demonstrated the existence of 
the strong coupling regime with an Anderson localized 
mode in two dimensions. The frequency splitting and 
its dependence on the parameters of the probe scatterer 
are described quantitatively using the coupled-mode the- 
oretical model, in which the parameters of the Ander- 
son localized mode are extracted from a spectrum of the 
LDOS. 

We shall show that an alternative formulation of the 
strong coupling criterion can be given, that is particu- 
larly relevant in the case of Anderson localization. Let 
us introduce the average linewidth of the electromagnetic 
modes 5uj and the average mode spacing Aoj. Normal- 
ized linewidths F^ = F;5 /Aa; and Tm = Em/5u) can be 
introduced, for the probe scatterer and for the localized 
mode M . F^ = 1 means that the bandwidth of the scat- 
terer covers on average only one mode of the disordered 
medium (the linewidth of the probe scatterer or emitter 
can be chosen or tuned to satisfy this condition). When 
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FIG. 2; (Color online) (a) Sketch of the system with a res- 
onant probe scatterer placed at the center, (b) Spectra of 
the dipole moment |ps(ij;)| of the probe scatterer for differ- 
ent values of the radiative linewidth Fg (from top to bottom 
rf = 10^-^ 25x10*^ s-^ SOxlO** s-^ 75x10® s''^; 10* s"^). 
(c) Frequency splitting in the spectrum of the dipole moment 
versus Ff. Solid line: Theoretical prediction by Eq. (|ll|l . 
Circles: Numerical simulations. 



the probe scatterer is resonant with localized mode M 
the strong coupling condition given by Eq. (|12p becomes 



1 Tm 
> — q , 



(15) 



where g = Sui/Auj is the normalized Thouless conduc- 
tance, a key concept in the theory of Anderson local- 
ization [23|. The localized regime corresponds to g < 1 
[this condition describes statistically a spectrum as that 
in Fig. Wic)]. The inequality shows that the smaller the 



conductance, the smaller the critical Purcell factor per- 
mitting to enter the strong coupling regime. This con- 
firms the idea that deeply localized modes in 2D or quasi- 
ID fis!, [l^ are particularly suitable to achieve strong 
coupling in the optical regime in condensed matter. For 
Tg ~ 1 and rj\/ ~ 1 (this condition is satisfied on average 
for the localized modes), the strong coupling condition 
takes the remarkable simple form Fp > g/2. This sim- 
ple relation directly connects the Purcell factor (a cen- 
tral quantity in cavity QED) and the Thouless conduc- 
tance (a statistical concept in transport theory). Finally 
let us remark that the inverse of the Thouless conduc- 
tance is statistically the analogue of the finesse of a stan- 
dard Fabry-Perot cavity that enters standard cavity QED 
analyses Q. 

In conclusion, we have demonstrated numerically the 
strong coupling regime between a resonant scatterer (or 
emitter) and an Anderson localized mode in two dimen- 
sions. The numerical results are in perfect agreement 
with a coupled-mode scattering theory in which the pa- 
rameters of the localized mode are extracted from a spec- 
trum of the LDOS. We have also shown that the strong 
coupling threshold can be expressed in terms of the Thou- 
less conductance and the Purcell factor. This analysis 
connects key concepts in transport theory and cavity 
QED, and might inspire the design of samples for the ob- 
servation of strong coupling to Anderson localized modes 
in photonics, or with other kind of waves. 
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